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We study the electronic properties of Mobius graphene strip with a zigzag edge. We show that 
such graphene strip behaves as a topological insulator with a gapped bulk and a robust metallic 
surface, which enjoys some features due to its nontrivial topology of the spatial configuration, 
such as the existence of edge states and the non-Abelian induced gauge field. We predict that 
the topological properties of the Mobius graphene strip can be experimentally displayed by the 
destructive interference in the transmission spectrum, and the robustness of edge states under certain 
perturbations. 
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I. INTRODUCTION 

Because of the unusual properties and potential ap- 
plications, topological insulators have recently been un- 
der great focus both experimentally [1, 2] and theoreti- 
cally 0, 0, H, 0, B, l8|, l9|l. The topological insulator sys- 
tem belongs to a novel category, possessing an insulating 
bulk, with a gap in the energy spectrum of propagat- 
ing electrons, whereas its surface is metallic. The edge 
states promise such metallic feature and describe the elec- 
trons localized on the surface with energy levels lying 
just within the gap of the bulk [s*]. For many systems in 
this category, the surfaces are no longer conductive when 
some perturbations are applied. However, for some sys- 
tem with nontrivial topology, the edge states are robust 
under perturbations. Such systems with robust metallic 
surfaces are referred to as topological insulators 01 • 

It is natural to imagine that those topological features 
of electrons can be realized through nontrivial topology 
in the configuration space of the system considered. A 
most recent illustration is the tight binding model for 
electrons hopping on a Mobius ladder [10]. In this inves- 
tigation, observable effects of Mobius boundary condition 
were found for a finite lattice. Destructive interference 
emerges from the transmission spectrum to display the 
typical topological feature. It was proved that such de- 
structive interference can be explained in terms of a non- 
Abelian gauge field induced by the nontrivial topology. 
However, such a novel nanostructure can not be regarded 
as a topological insulator, because this quasi- ID system 
has no edge states. 

In this paper, we study the electronic properties of 
a Mobius graphene strip. We show that with a 2D 
nontrivial topological structure, the Mobius graphene 
strip fTH, [ii, H Q E [iB'l with zigzag edges (see 
Fig-Hfb)) behaves as a typical topological insulator since 
it possesses robust edge states. It is noticed that there 
is no edge state in such Mobius strip with an armchair 
edge, thus no such nontrivial topological properties ap- 
pear. Our investigation will focus on the situation with 
a zigzag edge. Through analytic approaches, we first 
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FIG. 1: (color online) (a) A zigzag graphene strip is schemat- 
ically illustrated in a 2D version with open boundary con- 
dition along y direction. The white (black) dots represent 
sublattice A (B). The 3D version of the graphene strip with 
Mobius boundary condition along the x direction ('0^(5) (r) = 
'0b(a)(r+ivex), where r = (x,y) and? = (x, —y)) and periodic 
boundary condition {ipa(b){'^) = '0a(b) (r + Lex)) are plotted in 
(b) and (c), respectively. 



predict that the robustness of the edge states is main- 
tained under perturbation with uniform electric field in 
our discussion. And we compare the Mobius graphene 
strip with a generic one whose edges turn to be insulators 
in the same electric field. Besides, destructive interfer- 
ence in the transmission spectrum is found in the Mobius 
graphene strip, which is caused by the non-Abelian gauge 
field. Based on this observation, we propose a possi- 
ble approach of quantum manipulation on the transport 
properties of the Mobius graphene strip through a mag- 
netic flux. 



This paper is arranged as follows. In Sec. [Til we 
give an analytical description of the edge states for the 
generic and Mobius graphene strip. Here the tight bind- 
ing Hamiltonian is used to model the graphene strip with 
finite width. In Sec. lIIH a uniform electric field is applied 
to the generic and Mobius graphene strip as a perturba- 
tion. We compare the energy bands of the edge states in 
the generic and Mobius graphene strips under the uni- 
form electric field. Robustness of the edge states in the 
Mobius graphene strip under such perturbation is demon- 
strated. And we prove the candidacy of Mobius graphene 
strip as a topological insulator. Additionally, in Sec. HVl 
we discuss the topology-induced non-Abelian gauge field 
in the Mobius graphene strip and its observable effects in 
future experiments. Besides, we propose a possible quan- 
tum manipulation mechanism upon the Mobius graphene 
strip through a magnetic fiux. Finally, the conclusion is 
presented in Sec. [Vl 



II. EDGE STATES IN MOBIUS GRAPHENE 
STRIP 
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FIG. 2: (color online) On the left is the schematic illustration 
of the on-site electron occupation of an edge state with kx — 
1.27r/d, d— a/3/, in which the size of each circle represents the 
magnitude of the electron occupation on that lattice site. And 
the right figure displays the corresponding values of electron 
occupation on sublattice A (blue circles) and B (red circles) 
respectively. 



To describe the motion of electrons of graphene, we 
use the tight binding model [3, 



(1) 



where ?/^a(ri)('^6(ri)) annihilates an electron on site v\ of 
sublattice A(B), and the sum is taken over the nearest 
neighbor sites (z, j) with corresponding hopping constant 
J. For the zigzag graphene strip with finite width {N reg- 
ular hexagons in Fig.[T](a)) along y direction, the y com- 
ponent of the spacial vector r for A and B sublattices 
only take a finite number of values, e.g., y^ = [3{N — 
m) - l]l/2 for sublattice A, and yln^ = [3{N - m) - 2]//2 
for sublattice B. Here, / is the distance between nearest 
neighbors in the lattice, and m = 0, 1, ...,2A/' — 1. Usu- 
ally, it is assumed that the length of the graphene strip 
(M regular hexagons along x direction, as in Fig. Ufa)) 
is much larger than the width, namely, M ^ and 
periodic boundary condition is taken as 



V^a(6)(r) = V^a(6)(r + i^ex), 



(2) 



with L = >/3ML Thus, it is proper to perform Fourier 
transformation only for the x component of r upon the 
field operators of electrons as 

^a{h){k:,,y) = -^^V^«(5)(x,y)e-^^-^, (3) 



and kx can be considered as continuous since L is large 
enough. 

For a zigzag graphene strip, there exist eigenstates 
strongly localized on the edges of the strip with their 
energies lying exactly within the gap [H, [lO, [U, [22| . 



which are called edge states. A straightforward calcula- 
tion gives approximate edge states 



\^±(kx))=Dl^(kx)\vac), 



(4) 



which are defined by the corresponding annihilation op- 
erators 

De±{kx) = -^[Ae(kx) ± Be(kx)l (5) 
where the collective operators 

2N-1 



(6a) 



2N-1 



2N-l-m„ 



Mk.,yt^) (6b) 



m=0 



respectively represent parts of the edge states localized 
on the edges, as shown in Fig. O Here, 



Pk^ = -2cos 



V^kJ 



(7) 



represents the decay rate of the single electron on-site 
occupation probability along y direction. And 

is the normalization constant. \'^±{kx)) are the anti- 
bond and bond states with respect to the states 
Al{kx) \ vac) and Bl{kx) \vac). They correspond to the 
eigenenergy 

Ee±{kx)=± r'Y ' 
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FIG. 3: (color online) Schematic illustration of the Mobius 
graphene strip (a) before the unitary transformation W and 
(b) after the unitary transformation W. After the transforma- 
tion, the strip is divided into two independent pseudo strips 
(red t and blue | in (b)) with no interaction between each 
other. Besides, the sites on the pseudo edges y — ±1/2 obtain 
extra on-site potential, i.e., +J for the y = —1/2 sites and -J 
for the y = 1/2 sites. 



Actually, only when \pk^ \ < 1 are the edge states de- 
fined in Eq. ([5|) localized at upper and lower edges. The 
on-site electron occupation decays exponentially when 
heading into the bulk, which is plotted in Fig. O The 
requirement \pk^ \ < 1 means that kx can only be taken 
between the two neighboring Dirac points, namely, kx G 
(27r/3A/3/, 47r/3A/3/). No edge states exist beyond this 
region. 

It is pointed out here that only in the large N limit 
{N oo) are the above description for edge states 
(Eq. ([5])) accurate. However, for finite N, when kx is 
not in the vicinity of either Dirac point, namely 

- bfej) » 1, (10) 
the above description deviates merely negligibly from the 
I 



with 

,/.(ri,rj) = ^(ri-rj)-ex%i), (16) 

and 0{y) is a step function. Here the sum is taken 
over all the nearest neighbors except those whose bonds 
go across the ^ = line. 

From the above Hamiltonian we notice that the 
Mobius strip is divided into two separate generic pseudo 
strips, including y > strip (lower strip |) and y < 
strip (upper strip t) (see Fig. [3]). There is no cou- 
pling between the upper strip and the lower one, and 
the y = ±1/2 sites become "pseudo edges" after the 



accurate one. 

On the other hand, for a Mobius graphene strip with 
N regular hexagons along y direction and M (M N) 
along X direction (Fig. [T](b)), the Mobius boundary con- 
dition is explicitly written as 

'0a(6)(r) = V^5(a)(r + i^ex), (11) 

where r = (x^y) and r = (x, — ?/), and Gx is the unit vec- 
tor along X direction. To make Fourier transformation on 
X direction still available, a position-dependent unitary 
transformation W: ilJa,b i^a,(3 

Mr) =^IM - ^6(?)]e**^ (12a) 

Mr)=±[Mr)-Mr)W^'' (12b) 
for y > 0, and 

Mr) =^[Mr) + Mr)], (13a) 

Mr)=±[Mr)+M'^] (13b) 

for y < is necessarily used. It can be verified that after 
the transformation the new field operators ?/^Q!(^)(r) of 
the electrons satisfy the periodic boundary condition 

i^ai(3) (r) = i^a{(3) (i* + Le^)- (14) 

Then the Hamiltonian of the Mobius graphene strip be- 
comes H = Ho ± Hi, where 



(15a) 

(15b) 
I 

transformation. We should point out that these pseudo 
edges are not real in the spacial configuration of the 
Mobius strip. Actually, it follows from Eq. (|15ap and 
Eq. (|15bp that there is no coupling between V^q;(x, —1/2) 
and ijjf3{x,l/2) in H while there exist "on-site potential" 

on the pseudo edges ±J?/;j^(^^^j(x, =f //2)?/^q!(/3)(^7 tV^)- It 
has been proved that no edge states exist on these pseudo 
edges [23]. In fact, these pseudo edges appear because of 
the difference in unitary transformation between ^ > 
and y < and possess no topological features of real 
edges. 

As a consequence, there is merely one edge state for 



^0 = - J E ' {^i^OM^j) exp[^(/> (ri, rj)] + He.} - J^[V^l(x, -^-)M^, ^ + H.c], 



<i,j> 



Hi =J ^[V^i(x, -^)V^;3(^, b + H.c. + il;l{x,-^-)tlJa{x,-^-) - V^^(x, ^)V^/3(x, ^) 



4 



each X direction momentum kx respectively on both up- defined as 
per and lower strips, whose annihilation operators are 



J 



2N-1 



V^k^ m=N 



N-1 



(17a) 
(17b) 



with corresponding energies 

Ee^{kx)=J^^^^^. (18a) 

Eei{kx)=-Ee^{kx-^), (18b) 

respectively. 

Edge states in the graphene strip play an important 
role in the electron transport. For the monovalent model 
of a graphene strip, the Fermi energy is £^ = 0, around 
which edge states are highly degenerate. Besides, the 
energy band of the edge states is only half filled and the 
graphene strip is a zero-gap conductor. 
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III. A UNIFORM ELECTRIC FIELD APPLIED 
AS PERTURBATION 

In fact, the conducting behavior of the graphene strip 
(Mobius or generic) under perturbations (e.g., caused by 
a uniform electric field), is mainly determined by the 
properties of its edge states. For the generic graphene 
strip, a uniform electric field applied on y direction will 
cause a perturbation Hamiltonian 

He = -e^Et/4Vi(ri)V'a(ri) - 4{ri)M^i)], (19) 



where E is the electric field intensity. We assume here 
that the energy difference introduced by the electric field 
on opposite edges of the strip is much smaller than the 
hopping energy in graphene, namely, 3eENl <C J. Be- 
cause the energy difference between edge states and bulk 
states with the same kx is much larger than the pertur- 
bation, transitions between them, induced by the electric 
field, can be neglected. Thus we only focus on transitions 
between different edge states in the following discussion. 
Transition matrix elements between edge states are 

{^±{k,)\HE\^±{K))=0, (20a) 
{^±{K)\ He l^TiK)) =ekj{k:c, k',), (20b) 



FIG. 4: (color online) (a) Energy spectrum of a graphene 
strip (Mobius or generic), whose {N,M) = (10,100), with- 
out external electric field. Obviously the edge states between 
the two Dirac points (red line) are approximately degener- 
ate in the energy spectrum. A schematic illustration of the 
change in the energy spectrum is provided in (b) for a generic 
graphene strip and (c) for a Mobius graphene strip, which 
is caused by the applied uniform electric field on y direction 
whose intensity E satisfies eEl{3N — l)/2 = 0.15 J. Here the 
red solid line represents the energy level of the edge states 
without the electric field, while the green dashed line and the 
blue dot-dashed line are the changed energy bands of the edge 
states in the electric field. Here, as the dual vector of kx in 
the inverse Fourier transformation applied only on the field 
operators of the edge states, q is analogous to x. 



where \'^±{kx)) are edge states for the generic graphene 
strip defined in Eq. (|4]), and 



-eEl 




, (20c) 



with 5{kx^k'^) the Kronecker delta function. Thus the 
modified energies of the edge states are 



(21) 



5 



When is sufficiently large, 

eu.^e^^=-\eEl{m-l) 

is approximately independent of momentum kx. Equiva- 
lently, is approximately a constant in the range of 
where edge states exist. Besides, when the electric field 
is sufficiently strong but E <C J/3eNl still holds, the en- 
ergy Ee±{kx) of the original edge states approaches zero 



in comparison with |e/c^|, namely, \ek^ 
this sense, we approximately obtain 



i?:±(fc.)=±e 



(22) 



The above argument means that the originally highly 
degenerate energy level of edge states {Ee±{kx) — 0) is 
split into two separate energy levels by the uniform elec- 
tric field (see Fig. S^b)). This result could be simply 
interpreted by the different electric potential on the up- 
per and lower edges. Since a gap exists between the two 
energy levels, the zigzag edges of the graphene strip are 
no longer conductive. 

We notice that the above conclusion is only valid for 
the case with generic strips. For a Mobius graphene strip, 
due to its inherent twisted structure, the perturbation 
Hamiltonian introduced by the uniform electric field on 
y direction reads as 



H^E^^ = -eE^^4V^i(ri)V;a(ri)-V^J(fi)V^5(ri)]cos(^x,) 



eE 



2N-1 



E yL"M[V'i(fc.,yL"^)^/?(fc.,-yl"^) + ^f-c.] + [i^Uk, + ^,yL"))#(fc.,-t/L»)) + if.c.]}. (23) 



kx ,m=0 



L 



r 



Here we have assumed a uniform twist for the structure 
of the strip. It follows from Eq. ([23|) that the electric field 
couples the upper strip with the lower one, and induces 
transitions between states with different kx, making kx 
no longer conserved. For the same reason for a generic 
graphene strip, we still neglect the transition between 
edge states and bulk states. As a consequence, we only 
need to focus on the subspace spanned by all the edge 
states 



(24) 



The transition matrix elements between edge states on 
the same pseudo strip vanish, i.e.. 



^Ta)fe)|i^i^^|^Ta)(^^)) = o. 



(25a) 



while the nonzero matrix elements 



^,{k^J\ \^,{kx)) = {^,{kx)\ H^i"^ \^,{k',)) = ei''^\5{kx, fc;) + 5{kx + ^, fe;)] 



(25b) 



describe transitions between edge states on the upper 
and lower pseudo strips. Here e^^^ = e/e^/2, and we have 
neglected the difference between pk^ and Pfe^+27r/L- Still, 
we approximately obtain 



(M) _ (M) 

ti t n 

'^^ Vsi 



1 



eEl{3N-l). 



Since L is large, kx is approximately continuous, and a 
large number of kx that satisfy Mobius boundary condi- 
tion exist between the two Dirac points. The conclusion 
that ^eT/i(^a^) — is still valid for the Mobius graphene 
strip when kx is not in the vicinity of either Dirac point. 
Thus, the energies of the edge states in the presence of 



the uniform electric field are 



EL 



(26) 



where q is dual vector of kx in the inverse Fourier trans- 
formation. It could be concluded that the originally 
highly degenerate energy level Ef.y^{kx) = is broad- 
ened, by the uniform electric field, into an energy band 
with width 4e^^^^ (shown in Fig. [D^c)). A straightfor- 
ward explanation for the energy band broadening is that 
the electric potential on the edge of the Mobius strip 
varies along the x direction. 

It is important to point out that there exist no en- 
ergy gap for a Mobius graphene strip. The zigzag edge 
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FIG. 5: (color online) The density of state (DOS) 
e = is obtained by deducing the retarded Green's f 

tion G{e) = {e — H -\- ir])~'^ of the graphene strip 

T] = O.OIJ and (N, M) = (10, 100), and determining the spec- 
tral function A[e) — — 2ImG'(e). The black solid line, the 
green dashed line and the blue dash-dotted line represent the 
DOS of a graphene strip (Mobius or generic) without electric 
field, a generic graphene strip with uniform electric field, and 
a Mobius graphene strip with uniform electric field, respec- 
tively. 



is still conductive even in presence of an external electric 
field. A Mobius graphene strip with such typical feature 
is referred to as a topological insulator. Numerical calcu- 
lations for the density of states (DOS) near e = (Fig. [5]) 
obviously support the above theoretical predictions. The 
DOS-energy curve explicitly demonstrates the existence 
of edge states in the vicinity of e = when no external 
electric field is applied. In this case, for both the generic 
and Mobius graphene strips, there is only one peak cen- 
tered exactly at e = 0, indicating high degeneracy of 
edge states. When uniform electric field is applied on a 
generic graphene strip, the peak turns to two peaks cen- 
tered at ±e^^^p respectively, meaning that the highly 
degenerate energy level is split by the electric field. In 
this case, the DOS at e = approximately equals to zero, 
indicating the existence of a gap. Contrarily, with elec- 
tric field applied on a Mobius graphene strip, the DOS 
has apparently nonzero value for any e between ie^^^g^. 
This numerical result agrees well with our above analyt- 
ical prediction about the energy spectrum of generic and 
Mobius strips with and without external electric field, 
and indicates that the Mobius graphene strip is a topo- 
logical insulator, but a generic strip is not. 



IV. NON-ABELIAN GAUGE FIELD AND 
OBSERVABLE EFFECTS IN TRANSMISSION 

In this section, we study the transportation proper- 
ties of electrons in the Mobius graphene strip to demon- 




FIG. 6: (color online) Schematic illustration of a Mobius 
graphene ring with two leads connected to its opposite sides 
and an external magnetic flux thread its center. 



strate the topological effects similar to that for Mobius 
ladders [10]. By comparing the Hamiltonian Eq. (|15ap of 



a Mobius graphene strip with that of a generic graphene 
strip, it is recognized that, in the Mobius graphene strip, 
there is a phase shift on the hopping constant for the 
nearest neighbors {i^j) whose relative vectors ri — rj have 
nonzero x component. Similar to the case of the Mobius 
ladder, this phase shift in Eq. ([16]) can be described in 
terms of a non-Abelian gauge field 

A(r) = (A,,A,) = (^,0)%) 

in the continuous limit. The presence of this gauge field 
changes the canonical momentum p into p + eA(r). The 
gauge field only exists in the lower strip, to result in 
an effective magnetic flux $ = ch/2e along the positive 
y direction thread the center of the lower strip, which 
is bent into a ring. Thus when an electron travels one 
round along the lower strip, such a gauge field may bring 
about a phase shift tt. This effect could be experimen- 
tally testable when leads are connected on opposite sides 
of the strip (see Fig. [6]). Then, there would be no trans- 
mission between two leads through the lower strip, due 
to the destructive interference between electrons passing 
through two possible paths from one lead to the other. 
Such destructive interference is analogous to that in the 
usual Aharonov-Bohm effect [H, [2^1. However, this phe- 
nomenon in our system is totally induced by non-trivial 
topology, instead of a real magnetic flux. 

To demonstrate the non-Abelian nature of the gauge 
field more clearly, we only consider the special graphene 
strip with only one regular hexagon along y direction 
{N = 1), which is simply an aromatic hydrocarbon chain 
(see Fig. [6]). There are 4 possible values of y in the lattice 
of this strip, y = ±//2 and y = ±1 respectively. The 
spinor representation of electrons in a regular hexagon 
lattice is 



^(x) 



2 



V i^b{x ■ 



(27) 



-I)/ 



7 



where we have neglected differences in the x coordinates 
between field operators and the y coordinates have been 
omitted. Performing the unitary transformation W on 
^(x) , we have 



^{x) = W^{x) = — 
v2 



{M^)-^',{x))e^~ 



(28) 



It can be verified that ^(x) is a periodic function since 
^{x ^ L) = ^{x). Then the Hamiltonian of the Mobius 
graphene strip becomes 

H = -j"^^'^{x)U^{x)^[^'^{x)T^{x^V3l)^H.c.] 

X 

= -jJ2^Hkx)[U + {Te'^''''^ + H.c.)]^{ka:), (29) 
where 
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FIG. 7: (color online) (a) and (b) show the transmission spec- 
tra of a generic graphene strip and a Mobius graphene strip, 
respectively, both of whose parameters (N,M) = (1,40). 
Comparing (b) with (a), we find that several peaks in the 
transmission spectrum of the generic graphene strip are miss- 
ing in that of the Mobius graphene strip, due to the existence 
of a non-Abelian gauge field in the latter one. 
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and 



In the continuous limit, the above Hamiltonian be- 
comes 

H = [ d^^\^){v[R{-i^)-Ao]^Ro}^{^), (30) 
Jo 



dip 



where (p = 2itx/L^ v = 27rV^/J/L, R = diag[— a^y, a^y] 
{ay is the Pauli matrix), Aq =diag[cr^/2, 0], and Rq = J 
diag[0, 1,-1,0]. We would like to point out that in the 
continuous limit, we have expanded the original Hamil- 
tonian in the kx representation around kx = 7r/\^l. Be- 
cause Aq does not commute with Rq , Aq is regarded as a 
non-Abelian gauge field. Actually, the observable effect 
of the non-Abelian gauge field can be explicitly demon- 
strated when the strip has small width, since in this case, 
a graphene strip may have distinct energy bands instead 
of strongly overlapped ones. 

To obtain the transmission spectrum, we calculate the 



self energy [24 | 
^L{kx^ kx) = 



J' 



/ 1 





\ 1 



2 



/ 1 \ 





V 1 y 



(31a) 



(31b) 



in the matrix representation (actually the whole matrix 
is a 4M X 4M one made up of above 4x4 blocks 
with respective (/c^, kx))^ and J' is the hopping energy in 
the leads {J' = 1.5 J in our numerical calculation). The 
above self energy results from the connection of two leads 
to the strip. Then we obtain the effective retarded Green 
function of the graphene strip 



G{Eik)) = 



1 



E{k) - H + ir]-EL -^i 



(32) 



where E{k) = —2 J' cos{V^kl) is the energy of electrons 
injected through the left lead, and 7^ = 0+. After we de- 
termine the level broadening matrix Tl^r = —21m.J^L,Rj 
the transmission coefficient 

T{E{k)) = Tr[TRG{E{k))TLG\E{k))] (33) 

is obtained in a straightforward way. The transmis- 
sion spectrum of a generic graphene strip and a Mobius 
graphene strip (both width are = 1) are displayed 
in Fig. [71 (a) and (b), respectively. It is illustrated by 
Fig-E^b) that incident electrons with energy E G [— J, 0] 
are totally refiected by the Mobius graphene strip. This 
numerical result just confirm our heuristic prediction. 
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FIG. 8: (color online) (a) and (b) show the transmis- 
sion spectra of a Mo bins graphene ring, whose parameters 
{N, M) = (1, 40), with an external magnetic flux = n$o /2 
thread its center, where n is an integer, (a) corresponds to 
the n =even case, and (b) corresponds to the n =odd case, 
between which there are obvious differences in the position of 
the peaks in the transmission spectrum. 




(a) (b) 

FIG. 9: (color online) Schematic illustration of quantum 
transmission in a Mobius graphene ring with two leads con- 
nected to its opposite sides and an external magnetic flux 
$j = n$o/2 thread its center. Here, the up and down ar- 
rows represent the edge states in the upper and lower pseudo 
strips of the Mobius strip (see Fig. [3]). (a) When n =even, 
the only current carriers are electrons occupying upper edge 
states above the Fermi level (red dots), (b) When n =odd, 
the only current carriers are holes in the energy band of lower 
edge states below the Fermi level (light yellow dots). 



Since the non-Abelian gauge field exists in the lower 
strip instead of the upper one, there is a possible way 
to manipulate the transmission properties of the Mobius 
graphene strip through an external magnetic flux. As 
a magnetic flux is applied thread the center of the 
Mobius strip (which has the shape of a ring), a mag- 
netic vector potential appears on both the upper and 
lower pseudo strips. When the magnetic flux is half 
of integer magnetic flux quanta, i.e., = n<l>o/2 with 
$0 = ch/e (n = 1, 2, ...), the total effective magnetic flux 
in the lower strip becomes (n + l)<l>o/2 and in the upper 
strip it becomes n<l>o/2. Consequently, when n is even, 
quantum transmission is still suppressed in the lower 
strip, due to the destructive interference of electrons 
passing through two possible paths along the Mobius 
ring. Without such destructive interference, quantum 
transmission is allowed in the upper strip in this case. 
However, when n is odd, quantum transmission is al- 
lowed in the lower strip and suppressed in the upper one. 
Based on the above discussion, changes in the positions 
of the peaks in the transmission spectrum are expected 
to be experimentally observed when the external mag- 
netic flux changes. Numerical calculations, illustrated 
in Fig. [HI clearly verify our above heuristic discussions. 
Therefore, this magnetic-flux-based operation obviously 
implements quantum manipulation for electron transport 
in this Mobius nanostructure. 

Next, we further consider the quantum transport of 
electrons in the low energy excitation regime. Actually, 
these electrons locate in the energy band of the edge 
states. The energies of edge states in both the lower 
and upper strips (see Fig. [3]) are close to the Fermi level 
of the system. Here, the energies of the lower edge states 
are below the Fermi level, thus these edge states are all 
occupied at zero temperature. Contrarily, the upper edge 



states are not occupied at zero temperature, since their 
eigenenergies are above the Fermi level. 

When there are an integer number of flux quanta 
thread the Mobius ring (n=even in the above discussion), 
or there exists no external magnetic flux thread the ring, 
the fact that quantum transmission is forbidden in the 
lower strip means that holes below the Fermi level cannot 
be current carriers. Then the only current carriers in the 
Mobius ring between two leads are electrons occupying 
the upper edge states (see Fig. (Ua)). Oppositely, when 
the number of flux quanta thread the ring is half-integer 
(n=odd in the above discussion), quantum transmission 
is forbidden in the upper strip. In this case, holes of the 
lower strip right below the Fermi level, instead of elec- 
trons occupying the upper edge states, become current 
carriers between the two leads (see Fig.lSJb)). Moreover, 
if the electron-electron and electron-phonon interactions 
could not be ignored in the Mobius graphene strip, the 
transmission rate of electrons and holes would be differ- 
ent significantly. Thus, accompanying the switch of cur- 
rent carriers between electrons and holes, transmission 
rate in the Mobius graphene strip may change as well. 

Finally, we remark on the effect of electric field on the 
transmission properties of the Mobius graphene strip. 
The uniform y direction electric field applied to the 
Mobius graphene strip can induce strong coupling be- 
tween electrons in the lower and upper strips. In this 
case, both of the strips can contribute significantly to 
quantum transmission. Therefore, the above forbidden 
transmission in one of the two pseudo strips no longer 
emerges in this case. 
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V. CONCLUSION 

Oriented by physical realizations of topological insu- 
lators and topological quantum devices, we theoretically 
studied the electronic properties of the Mobius graphene 
strip, which is an exotic 2D electron system with a topo- 
logically non-trivial edge. Various properties of the edge 
states were investigated through the tight binding model 
in this paper. We also studied the robustness of edge 
states in the Mobius graphene strip under perturbations, 
such as that caused by a uniform electric field. Analytical 
results about the exotic natures of such electron system 
are obtained for the first time and then confirmed by the 
following numerical calculations. Moreover, the physi- 
cal effects of the non-Abelian induced gauge structure in 
the Mobius graphene strip were studied by considering 
the transmission spectrum of the Mobius graphene strip 
with two leads connected to its opposite sides. Based on 
such non-Abelian induced gauge field discovered, which is 



similar to that for the Mobius molecular devices studied 
in Ref. [10], we proposed a possible manipulation mech- 
anism for the Mobius graphene strip through a magnetic 
fiux. In fact, the robust edge of the Mobius graphene 
strip and the non-Abelian induced gauge field are two 
characteristic demonstrations of the nontrivial topology 
of the Mobius graphene strip. The above characters make 
the Mobius graphene strip as a candidate for topological 
insulators, which may benefit future applications in quan- 
tum coherent devices, and quantum information process- 
ing. 
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